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PARTICLE  BEAM  CHEMISTRY  STUDIES 


T  INTRODUCTION 

The  successful  propagation  of  a  pulsed  particle  beam  depends  in 
part  on  the  conditions  in  the  air  channel  which  is  formed.  In  particular, 
*  the  conductivity  is  a  critical  parameter,  and  this  is  determined  primarily 

by  the  electron  density.  This  report  describes  a  short  study  which  con¬ 
sisted  of  several  chemistry  calculations  in  the  channel  following  an 

intense  2  mm  radius  beam  and  one  calculation  of  the  chemistry  in  the  beam. 

I 

CHANNEL  INITIAL  CONDITIONS 

The  temperatures  and  species  concentrations  at  the  end  of  the 
beam  pulse  where  taken  from  a  calculation  made  by  SAI 1  for  a  10  kA,  2  mm 
radius  pulse.  Although  the  pulse  length  was  10-8  sec,  the  SAI  calculation 
was  carried  out  to  2xl0"8  sec  because  of  considerable  ohmic  heating  which 
occurs  after  the  high  energy  particle  beam  passes.  The  conditions  at 

I 

2xlo~8  sec  were  used  as  initial  conditions  in  our  calculation.  The  heavy 
particle  temperature  was  1810K,  the  electron  temperature  was  6380K  and  the 
nitrogen  vibrational  temperature  was  5642K.  The  species  concentrations 
are  shown  in  Table  1. 

The  case  chosen  by  SAI  was  a  particularly  intense  beam  and 
involves  the  deposition  of  a  considerable  amount  of  energy.  The  total 
energy  deposited,  including  ohmic  heating,  is  about  3.5  x  107  ergs/cm3. 
About  thirty  percent  of  this  energy  is  contained  in  N2  vibration  and 
thirty  percent  in  dissociation  and  electronic  excitation.  This  energy  can 
be  released  by  chemical  reactions  and  excited  state  quenching  during  the 
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Table  1.  Species  concentrations  at  2  x  10“ 8  sec. 
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*  Number  In  parenthesis  Is  the  power  of  ten. 


evolution  of  the  channel  and  can  significantly  Impact  the  gas  and  electron 
temperatures. 
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CHANNEL  EXPANSION 


The  heating  of  the  gas  will  cause  an  expansion  of  the  channel. 
The  initial  heating  of  the  gas  to  1810K  will  cause  an  expansion  (assuming 
adiabatic  expansion,  vibrational  freezing  and  no  further  energy  release) 
of  about  a  factor  of  four  and  a  concurrent  cooling  of  the  gas  to  about 
1070K.  This  expansion  takes  several  microseconds  to  occur  but  since  we 
will  be  interested  in  tens  of  microseconds  we  have  assumed  the  expansion 
takes  place  instantaneously.  We  therefore  reduced  all  the  species  concen¬ 
trations  in  Table  1  by  a  factor  of  four  and  used  these  values  as  our 
starting  concentrations.  In  addition  we  made  an  estimate  of  the  N20  which 
would  have  formed  during  the  energy  deposition  and  added  this  to  the 
initial  conditions.  This  was  done  because  N20,  at  these  temperatures,  is 
an  excellent  attacher  of  electrons.  The  value  of  N20  we  assumed  after 
expansion  was  3.6  x  1015  cm"3. 

CHEMISTRY  CALCULATIONS 

The  chemical  evolution  of  the  channel  was  calculated  using  our 
SEACHEM  code.  This  code  integrates  the  coupled  differential  chemical 
equations  and  contains  about  70  species  and  about  600  reactions.  Since 
the  chemical  development  will  involve  the  release  of  energy  stored  in 
dissociation  and  vibration  during  the  energy  deposition  period,  the  code 
keeps  track  of  the  energetics  and  allows  for  heating  of  the  gas.  It  also 
allows  for  instantaneous  adiabatic  expansion  of  the  gas  as  heating  occurs. 
This  density  decrease  amounted  to  factors  of  between  1.5  and  2  in  our 
calculations,  depending  on  the  assumptions. 

We  will  not  list  all  the  reactions  in  our  code.  As  it  turned 
out  the  high  temperatures  in  the  channel  simplify  the  chemistry  and  only 
several  reactions  or  processes  are  critical  to  the  electron  density. 


When  the  ionization  source  turns  off  at  10~8  seconds  the  domi¬ 
nant  ion  is  02+.  Because  of  charge  exchange  reactions  N0+  quickly  becomes 
the  dominant  ion  thereafter  and  by  10"7  seconds  is  the  only  significant 
positive  ion.  The  electron-ion  recombination  rate  of  N0+  is  therefore  an 
important  rate  constant.  We  have  used  a  value  of 

kj  *  4*10”7  cm3/sec  ,  (1) 

^e 

where  Te  is  the  electron  temperature.  There  is,  however,  considerable 
uncertainty  in  the  dependence  on  Te.  Experimental  values  range  from 
Te"*4  to  Tg'1.  For  an  electron  temperature  of  about  0.5  eV  or  6000K, 
this  amounts  to  a  difference  in  the  rate  constant  of  a  factor  of  six.  In 
a  recombination  dominated  decaying  plasma  this  amounts  to  a  factor  of  six 
difference  in  the  electron  density. 

Electron  attachment  is  usually  an  important  parameter  in  deter¬ 
mining  the  electron  density.  The  main  attachment  mechanism  is  usually 
three  body  attachment  to  02,  i.e., 

02  +  e  +  M  ♦  02~  +  M  (2) 

The  production  of  N20,  however,  in  large  amounts  makes  the  following  reac¬ 
tion  potentially  important, 

e  +  N20  -►  NO  +  0"  .  (3) 

Although  this  reaction  is  endothermic  by  about  0.2  eV,  the  presence  of 
both  hot  electrons  and  vibrational ly  excited  N20  can  make  the  rate  con¬ 
stant  for  reaction  3  significant.  Using  the  experimental  data2  and  an 
electron  temperature  of  several  tenths  of  an  eV  we  have  obtained  the  rate 
constant 
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cm3/sec. 


(4) 
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k3  =  10-8  f_L)  e 

300 


-3200/T 


where  T  is  the  N20  temperature. 


Detachment  will  also  be  important  and  the  primary  detachment 


mechanisms  are 

02”  +  M  -  02  +  e  +  M  (4) 
02-  +  0  -  03  +  e  (5) 
O'  +  0  +  02  +  e  (6) 
0"  +  02( lA)  03  +  e  (7) 
02-  +  02(1a)  ♦  02  +  02  +  e  (8) 


Reaction  4  is  important  because  of  the  high  temperatures  involved  in  the 
channel . 

Ion-ion  recombination  is  usually  dominated  by  the  three-body 
process  at  sea  level  density  and  temperature.  8ecause  of  the  heating  and 
expansion  of  the  channel,  however,  the  two-body  process  must  be  included. 
The  effective  two  body  rate  constant  used  in  our  calculations  was 

kg  s  6xl0’8  (300)°  +  10-25  (3Q0')2*5  (m)  cm3/sec  (9) 

where  T  is  the  gas  temperature  and  (M)  is  the  total  particle  density  in 
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Finally,  because  of  the  energy  stored  in  N2  vibration,  the 
quenching  of  N2  vibration  is  an  important  parameter.  Atomic  oxygen  seems 
to  be  the  most  efficient  quencher  of  N2  vibration  of  all  the  species 
involved  in  our  calculations.  The  measured  rate  constants  between  1000K 
and  2000K  can  be  fit  by  the  expression3 

k10  =  7.8x10-15(_L)1,7  cm3/sec  (10) 

300 

A  question  which  arises,  however,  is  how  efficient  is  the  quenching  by 
oxygen  of  the  high  vibrational  states  of  N2.  If  the  quenching  of  the  high 
vibrational  states  can  also  be  given  by  the  above  expression  then  the  N2 
will  lose  a  significant  portion  of  it's  energy  to  translation  in  several 
tens  of  microseconds.  If  not  it  will  take  significantly  longer. 

RESULTS 

Four  calculation  were  made.  In  two  of  the  calculations  the  N2 
vibration  was  assumed  to  be  quenched  slowly  while  in  the  other  two  the 
quenching  was  assumed  to  be  rapid.  Since  the  interchange  of  energy  between 
electrons  and  nitrogen  vibration  is  very  rapid4  we  assumed  that  the  elec¬ 
tron  temperature  was  equal  to  the  N2  vibrational  temperature.  Thus  the 
quenching  of  N2  controls  the  electron  temperature.  In  each  pair  of  calcu¬ 
lations  one  calculation  assumed  no  water  vapor  and  the  other  assumed  a  3 
percent  H20  concentration.  The  H20  was  added  to  see  if  this  made  a  sig¬ 
nificant  difference  in  the  electron  density  due  to  clustering  of  the 
positive  ions  and  the  formation  of  ions  such  as  H20+  and  H30+.  In  addi¬ 
tion,  the  inclusion  of  H20  will  yield  species  such  as  H,  OH  and  H02. 
These  species  can  greatly  effect  the  oxygen  chemistry  and  therefore  impact 
on  the  nitrogen  vibrational  quenching. 


> 


The  electron  densities  for  these  four  cases  as  a  function  of 
time  (in  microseconds)  are  given  in  Table  2.  Cases  A  and  B  are  for  mini¬ 
mum  N2  vibrational  quenching.  Case  A  is  for  no  H20  and  case  B  is  for  3 
percent  H20.  Cases  C  and  D  are  for  maximun  N2  vibrational  quenching,  with 
case  C  having  no  H20  and  case  D,  3  percent  H20. 


Table  2.  Electron  densities. 


t(ysec) 

Case  A 

Case  B 

Case  C 

Case  0 

10 

2.1  (12)* 

2.1  (12) 

1.8  (12) 

1.9  (12) 

20 

1.1  (12) 

1.1  (12) 

1.1  (12) 

1.1  (12) 

40 

5.9  (11) 

6.1  (12) 

6.4  (11) 

6.3  (11) 

70 

3.8  (11) 

3.8  (11) 

3.5  (11) 

3.1  (11) 

100 

2.9  (11) 

2.7  (11) 

2.1  (11) 

1.9  (ID 

*  Number  in  parenthesis  is  power  of  ten. 


» 
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We  see  that  there  is  very  little  variation  in  the  electron 
density  among  the  different  cases  except  at  late  times  and  then  only  by, 
at  most,  50  percent.  The  addition  of  water  vapor  makes  no  significant 
difference  in  the  type  of  ion  which  is  present.  Because  of  the  high 
temperatures  involved  there  is  no  clustering  and  N0+  remains  the  dominant 
ion.  The  addition  of  H20  does  make  a  significant  difference  in  the  recom¬ 
bination  of  atomic  oxygen.  The  atomic  oxygen  disappears  much  more  rapidly 
and  this  does  decrease  the  N2  quenching.  However,  part  of  this  decrease 
is  offset  by  quenching  by  H  and  H2.  Also  the  more  rapid  recombination  of 
oxygen  adds  to  the  heating  of  the  gas  and  this  offsets  the  reduced  gas 
heating  caused  by  the  reduced  N2  vibrational  quenching.  A  more  accurate 
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treatment  of  the  interaction  between  the  electron  temperature  and  N2 
vibration  would  probably  make  the  differences  between  the  water  and  no 
water  cases  larger  but  they  would  still  not  be  significant. 

The  presence  o*  N20  does  increase  the  attachment  significantly 
at  early  times.  At  late  times  the  high  temperature  has  caused  much  of  the 
N20  to  thermally  dissociate.  Also  at  late  times  the  attachment,  which  is 
still  large  even  without  the  N20,  is  offset  by  an  equally  large  detachment 
caused  in  large  part  by  the  high  temperature.  Thus  the  electron  density 
decay  is  electron-ion  recomination  dominated.  The  electron  density  is 
essentially  determined  by  the  dissociative  recombination  rate  constant  for 
N0+. 

CALCULATION  WITHIN  THE  BEAM 

The  addition  of  water  vapor  in  the  previous  calculations  had  no 
effect  in  changing  the  fact  that  N0+  was  the  dominant  ion.  However,  the 
SAI  calculation  within  the  beam  had  no  water  vapor.  The  ions  N2+,  0+  and 
N+  can  charge  exchange  with  water  to  form  H20+  and  this  ion  can  react  with 
water  to  form  H30+.  We  thought  that  the  H30+,  which  has  a  larger  recom¬ 
bination  rate  than  N0+,  formed  within  the  beam  might  persist  and  contri¬ 
bute  to  the  electron  decay  after  the  beam  had  passed.  We  therefore  made  a 
calculation  within  the  beam.  The  ion-production  rate  specified  by  SAI  was 
used  (about  4  x  1017  ion-pairs/cm3/sec  for  10~8  seconds)  as  was  their 
results  for  the  gas  and  electron  temperatures.  The  calculation  showed 
that  H20+  and  H30+  were  indeed  significant  within  the  beam  but  that  these 
species  decayed  rapidly  after  the  ionization  source  is  turned  off.  By 
2  x  10“8  seconds,  the  recombination  rate  was  totally  dominated  by  N0+  and 
02+. 

This  calculation,  however,  did  show  significant  differences  in 
the  concentrations  of  a  nunber  of  other  species  at  2xl0*8  seconds  between 


the  SAI  calculation  and  ours.  Our  atomic  oxygen  concentration  is  about 
half  that  obtained  by  SAI  and  this  is  due  to  H  and  OH  chemistry.  Our 
02( *a)  concentration  was  only  40  percent  the  corresponding  SAI  value  and 
this  was  probably  due  to  our  inclusion  of  quenching  by  electrons  and 
perhaps  a  lower  value  of  the  02( xa)  formed  per  initial  ion-pair.  These 
reduced  values  could  impact  on  later  detachment  and  N2  vibrational  quench¬ 
ing,  particularly  in  less  intense,  lower  temperature  beam  cases. 

Our  N(2D)  concentration  was  about  an  order  of  magnitude  smaller 
than  the  SAI  value  and  this  is  undoubtedly  due  to  the  quenching  of  this 
state  by  H20.  This  quenching  probably  causes  the  species  NH  to  be  formed 
and  this  causes  a  whole  new  chemistry  problem  to  be  introduced  since  the 
chemistry  of  NH  is  largely  unknown. 

Our  value  of  ozone  was  considerably  larger  than  that  in  the  SAI 
calculation.  We  obtained  an  03  concentration  of  8.7xl014  cm"3,  which  was 
about  400  times  the  SAI  value.  This  large  formation  was  due  to  the 
detachment  reactions 


02"  +  0  +  03  +  e 

(11) 

0"  +  02( xa)  +  03  +  e 

(12) 

This  may  be  somewhat  unrealistic  since  the  reverse  of  these  reactions  were 
not  included  even  though  they  are  only  about  0.6  eV  endothermic.  The 
rate  constants,  however,  cannot  be  obtained  from  detailed  balance  since 
these  reactions  almost  certainly  involve  03  in  vibrational ly  excited 
states.  The  dependence  of  the  rate  constants  on  the  vibrational  states  is 
unknown  and  must  be  determined  experimentally. 

Finally  the  electron  density  we  obtained  at  2xl0-8  seconds  was 
8.3  *  10 14  cm"3,  a  factor  of  about  four  less  than  that  in  the  SAI 
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calculation.  This  is  due  to  the  attachment  to  N20  formed  in  our 
calculation.  As  we  have  already  pointed  out,  however,  this  N20  will  be 
destroyed  later  by  the  high  temperature. 

CONCLUSIONS 

The  main  conclusion  to  be  drawn  from  our  calculations  is  that 
the  chemistry  of  very  intense  beams  is  fairly  simple  due  to  the  high 
temperatures  involved.  The  electron  concentration  is  primarily  determined 
by  the  electron-ion  recombination  rate  constant  for  N0+.  Unfortunately 
this  rate  constant  is  highly  uncertain  at  high  temperatures.  We  used  a 
Te  *1  dependence.  Use  of  a  Te"0,lt  dependence  would  probably  reduce  the 
electron  density  to  below  10 11  cm3/sec  after  40  microseconds.  This  could 
have  serious  consequences  for  the  stability  of  the  second  pulse  trying  to 
propagate  in  the  channel  at  those  times.5  Another  uncertainty  resides  in 
the  electron  temperature  at  microsend  times  due  to  the  uncertainty  in¬ 
volved  in  the  quenching  of  N2  vibration  and  its  coupling  to  Te.  This 
should  be  examined  more  carefully. 

There  is  considerable  interest  in  less  intense  beams.5  In  these 
the  heating  will  be  much  less  severe  and  the  chemistry  in  the  channel  will 
be  more  complicated.  Cluster  ions  will  be  important  and  will  persist  out 
to  late  times.  Attachment  will  be  greater  due  to  N20  and  probably  03. 
This  should  be  investigated  in  more  detail.  Detachment  will  be  much 
less.  The  quenching  of  N2  vibration  will  be  reduced  and  this  will  keep 
the  electron  temperature  higher  for  a  longer  period  of  time.  Previous 
calculations  have  indicated  that  it  will  take  about  three  pulses  to  pro¬ 
duce  a  highly  ionized  channel  (ne>10n).  However,  in  these  calculations 
ohmic  heating  was  not  included,  Te  was  wrong,  N2  vibration  was  not 
included,  the  expansion  was  not  properly  handled  and  the  chemistry  scheme 
was  not  as  complete. 


REFERENCES 

1.  Keeley,  D. ,  and  R.  Johnson,  Private  Communication,  Stanford  Research 
Institute. 


2.  Chantry,  P.  J.,  OCP,  51,  3369,  1969. 

3.  McNeal,  R.  J.,  et  al.,  JGR,  79,  1527,  1974. 

4.  Ali,  A.  W.,  The  Physics  and  the  Chemistry  of  NRL  Master  Code  for  the 
Disturbed  E  and *T  Regions,  Nkl  Memorandum  Report  373z,  Naval  Research 
Laboratory,  l97£. 

5.  Yu,  S.,  Private  Communication,  Lawrence  Livermore  Laboratory. 


13 


